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On the basis of local nonequilibrium approach, the one-dimensional model of the solute diffusion 
during rapid solidification of the binary alloy in the semi-infinite volume is considered. Within the 
scope of the model it is supposed that mass transport is described by the telegrapher equation. 
The basic assumption concerns the behavior of the diffusion flux and the solute concentration at the 
interface. Under the condition that these quantities are given by the superposition of the exponential 
functions the solutions of the telegrapher equation determining the fiux and the solute distributions 
in the melt have been found. On the basis of these solutions different regimes of the solidification 
in the near surface region and the behavior of the partition coefficient have been investigated. The 
concentration profiles in the solid after complete solidification are analyzed depending on the model 
parameters. 



PACS numbers: 81.10.Aj, 05.70.Fh, 05.70.Ln, 81.30.Fb. 



I. INTRODUCTION 

In the present time the process of rapid solidifica- 
tion is a well established method for production of the 
metastable materials and, in particular, supersaturated 
solid solutions. The supersaturated metal structures can 
form during solidification of binary alloys due suppression 
of the processes of the solute segregation at the rapidly 
moving solid- liquid interface [lH3|. Quantitatively this 
effect can be characterized by the partition coefficient k 
defined as the ratio of the solute concentration in the 
growing solid to that in the liquid at the sharp interface. 
The phenomenon of "solute trapping" by the growing 
phase, implying, in particular, that the partition coeffi- 
cient deviates from its equilibrium value and increases 
towards unity at large growth rates, has been attracting 
considerable attention over of several decades both from 
experimental and theoretical points of view 4 23]. 

In the rapid solidification experiments very high veloc- 
ities of the phase interface can be reached so that the de- 
viations from local equilibrium in both the bulk of phases 
and at an interface become considerable 0, @, [3, IIS] ■ 
For theoretical description of solute trapping and related 
phenomena observed during rapid solidification a number 
models have been proposed [^-ly] in which, in particular, 
the deviation from local (chemical) equilibrium at solid- 
liquid interface is described by the partition coefficient 
k{V) depending on growth velocity V. According to the 
continuous growth model (CGM) of Aziz and Kaplan Q 
the velocity depending partition coefficient for dilute so- 
lutions is given in the form 



k{V) 



V/Vi 



DI 



1 + V/Vi 



DI 



(1) 



where Vdi is the atom diffusive speed at the interface. 

The equation H]) predicts that the value A: = 1, that 
is, the interface motion without solute partitioning, is 
only reached at V/Vdi ^ 1- Meantime there are a num- 
ber of the experimental works [l^, [25l - [28j in which it has 
been shown that the transition to complete solute trap- 
ping giving rase to diffusionless solidification occurs at 
substantially finite values of V/Vdi- This circumstance 
is automatically taken into account within the scope of 
the local nonequilibrium approach developed in the works 
p^ - [23| . At the high growth velocities the deviation from 
local equilibrium can be essential not only at the interface 
but in the bulk of the liquid phase as well. The extension 
of the expression ^ for this case introduced in [l7i] has 
the form 



k{V) = 



(1 - V^/V^)ke + V/Vdi 



1, 



V/Vi 



DI 



V <Vd 
V^Vd 



(2) 



The expression ([2]) takes naturally into account the fact 
that when the growth velocity V exceeds the velocity 
of the propagation of concentration disturbances in the 
liquid phase Vd the solute transfer in a melt is absent and 
the transition to the diffusionless solidification occurs at 
the finite velocity V — Vd ■ The extensions of expressions 
([T]) and ^ to the application to the concentrated alloys 
are given in [1, [2^ . 

It should be noted that the currently existing theoret- 
ical models consider the phase interface far from the sur- 
face of a melt that is in fact an infinite system. The pres- 
ence of the surface leads to the additional peculiarities. 
If the matter flux from the surface is absent, complete 
solute trapping by the growth phase must occur in the 
near surface layers. As a result, the partition coefficient 
must tend to unity when the interface approaches the 
free surface even with small velocities V and can vary in 
a complex enough way in the near surface region except 
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for the case oiV^Vo- 

The behavior of the key quantities, characterizing a 
sohdification process, at the moving interface specifies 
eventually the micro - and macrostructure of the result- 
ing solid and therefore is a subject of the central interest. 

Rapid solidification of dilute Ni(Zr) alloy after pulsed 
laser irradiation has been studied in the work of Arnold 
et. al. A simulation of the distribution of Zr after 

complete solidification using the CGM has shown that 
the simulated concentration profile reproduces the exper- 
imental data well enough in the deep regions. However, 
there is a significant discrepancy in data in the near sur- 
face region. It has been assumed in that an initially 
planar interface becomes unstable in the near surface re- 
gion and degenerates into a cellular structure although 
the authors have not been able to verify this experimen- 
tally. On the other hand, a one-dimensional simulation 
using a planar interface and the partition coefficient in- 
creasing with time in near surface region as 



kit) 



(3) 



where r is a fitting parameter, has given good enough 
agreement with the data in near surface too. 

The relationship ^ enables one to make some assump- 
tion about the behavior of the solute concentration at the 
interface. Indeed, according to the relationship ((3]) the 
dependence of the partition coefhcient on the position of 
the interface x = —Vt, moving with a constant velocity 
in parallel to the free surface (fixed dX x — Q) has the 
form 



k{x) 



zjl 



(4) 



with I = tV . Because k is defined as the ratio of the so- 
lute concentrations in different phases, the relationship 
(j4|) means in the simplest case that the solute concen- 
trations taken at the interface in each phase can vary 
exponentially as well but with different exponents. In a 
more general case these quantities can be defined as the 
superposition of the terms of the form ^ . 

In the present work the one-dimensional model of the 
solute diffusion during rapid solidification of the dilute bi- 
nary alloy in the semi-infinite volume is considered. The 
interface velocity is presumed to be high enough to pro- 
vide stability of the planar front of the growth of the solid 
phase [2l|, [2^ . It is also presumed that the motion of the 
interface is so fast that local equilibrium in the bulk of 
the liquid phase is absent and the solute transport oc- 
curs under local nonequilibrium conditions 'l6M23l]. Such 
conditions are achieved by high undercooling of the melt 
or during recrystallization after laser irradiation of the 
solid. The main purpose of this work is the development 
of the model determining the behavior of the solute con- 
centration and the diffusion flux both at the fast moving 
interface and in the bulk of the phases and the determi- 
nation of the inhomogeneous partition coefficient. 

The work is organized as follows. In Sec. [n]the general 
equations describing the local nonequilibrium transport 



are formulated and the model defining the behavior at 
the interface of the key quantities of interest is given. 
The solutions of the transport equations for semi-infinite 
volume have been obtained by and the Riemann method 
for the hyperbolic differential equations. These solutions 
define the diffusion flux in the liquid phase and the so- 
lute concentration in both phases. The inhomogeneous 
partition coefficient is derived. The discussion and the 
conclusion are respectively given in Sec. IIIII and Sec. IIVI 
The Riemann method in detail and its application to the 
present problem are given in Appendixes. 



II. THE MODEL 

As it has been noted above the interface velocity can 
reach high values. This takes place, for example, at the 
high undercooling of the melt or during solidification af- 
ter the laser irradiation of the surface of an alloy. When 
the interface velocity V is equal or more than the velocity 
of the propagation of the concentration disturbances in 
the liquid phase Vd, the diffusion fleld in the liquid can 
signiflcantly deviate from local equilibrium p^.Tl9j. In 
this case the solute diffusion flux is no longer defined by 
the classical Pick's law relating the diffusion flux to the 
gradient of a concentration, and it should be considered 
as an independent variable as well as the solute concen- 
tration. According to extended irreversible thermody- 
namics [m the simplest generalization of the Pick's law 
taking into account the relaxation to local equilibrium 
in the diffusion fleld is given by the Maxwell-Cattaneo 
equation which one writes down in the one-dimensional 
form as 



Jl 



T^T- = -Dl- 



dt 



dx 



(5) 



where index L relates to the liquid phase, Jl is the solute 
diffusion flux, Cl is the solute concentration, t is the time 
of relaxation of the diffusion flux to its local equilibrium 
value defined by the Pick's law and is the diffusion 
coefficient. 

Pq. ([5]) in combination with the conservation law 



dCL 
dt 



dJL 
dx 



gives rise to the hyperbolic transport equations 



d^CL 
^ dt^ 



dCL 
dt 

dt 



D, 



dx^ 



(6) 

(7) 
(8) 



The equation of the type (O and ^ is known as the 
telegrapher equation that combines the properties both 
of the wave equation and the diffusion one. At the times 
of the order t it predicts the finite velocity of the propa- 
gation of concentration disturbances Vd = (Dl/t)^^'^ in 
contrast to the diffusion equation for which Vd = oo at 
all time scales. 
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To describe the mass transport during the solidification 
process we consider a binary melt initially occupying half- 
space X ^ 0. The planar front of solidification forms in 
the infinitely removed region at i = — oo and isothermally 
moves with the constant average velocity V to the surface 
of the system, fixed at x — 0, along the trajectory x + 
Vt = in parallel to the free surface. At an arbitrary 
moment of time the region occupied by the melt is given 
by the inequality ^ x ^ ~Vt {t ^ 0). Therefore in the 
plane {x, t) the hquid phase occupies the region x + Vt ^ 
0,a;^0,i^O. At the interface representing the surface 
of a discontinuity the conservation law of mass holds [s^l 
that in accepted notations has the form 



should be fulfilled 



[Jl + VCl] 



x+Vt=(i 



[Js + VCs] 



x+Vt=Q 7 



(9) 



where index S relates to the solid. Taking into account 
a small mobility of the solute in the solid by compari- 
son with its mobility in the liquid phase one neglects, as 
usual, by diffusion in the solids and writes down Eq. ([9|) 
as 



V{Cl - Cs)\x+vt=o 



-J. 



L\x+Vt=0 



(10) 



Now we consider the diffusion flux in more detail. Intro- 
ducing dimcnsionlcss variables t/r, x/tVd in Eq. ([9|) one 
obtains 



9V 



dJ_ 

dt 



dx'^ 



(11) 



where the former notations {x, t) have been used for new 
variables and J/Vd is the dimensionless diffusion flux. 
The boundary condition (jlOp in the dimensionless form 
is written as 



Oi{Cl — Cs)\x+at=0 = — J|a:+at=0 



(12) 



where the dimensionless parameter a = V/Vd character- 
izes the extent of the deviation of the system from lo- 
cal equilibrium. In addition, at the surface the equality 



J(.Ti)U=o = 



(i ^ 0) , 



(13) 



expressing the condition of the absence of the flux 
through the surface. At last, the solution of Eq. (fTTI) 
is sought in the region X = x-|-ai<0, x^O, t^O 
occupied by the liquid phase while the solid occupies the 
region X ^ 0. 

It is physically apparent that at a ^ 1, that is V ^ Vd, 
the presence of the surface is not of considerable impor- 
tance. Interface moves with the velocity equal to or ex- 
ceeding the velocity of the propagation of the concentra- 
tion disturbances and the solute distribution in the liquid 
phase remains homogeneous. The solution of Eq. (fTTj) 
satisfying this condition and compatible with the equal- 
ity (ini) is [ii, m 

J = 0, Cl = Cs = const (a>l), (14) 
which corresponds to the complete solute trapping by the 
growth phase. 

Now we consider the case of a < 1. Suppose that at the 
moving interface residing in an arbitrary point x at the 
moment t — ~x/a the flux J and its the time derivative 
dJ/dt are known 



J{xt)\t 
dJ{xt) 



dt 



X I a. 

t — — X I OL 



^ jo (a;) , 



(15) 
(16) 



where the function jo (x) and ji (x) will be specified fur- 
ther. Eqs. (|15p and determinate the "initial" condi- 
tions that are given at the straight line x+at = defining 
the trajectory of the interface rather than at i = 0. 

If the functions jo{x) and ji{x) are known the solution 
of Eq. (fTTj) satisfying the conditions (fT5)) and (|16p in the 
region X ^ at a < 1 can be found by the Riemann 
method [33| (for details see Appendix A) and has the 
form 



a{x — t) 



X 



2(1 - a) 



X — t 
(p\a— — J cxp 

1 + a 



X 



2(1 + a) 



_ ip-t/2 
2 



dxi 



4a 



a{x -I- 1) 
1 — a 
a{x — t) 
l + a 



a{x + 1) 
1-a 



4(y{x-x,r~{t + x,/a)A 
^{x-xi)^ - (t + xi/a)2 



(17) 
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where 

^(a:)-jo(a:) (18) 

and Jq{x) is the Bessel function of zero order. At arbi- 
trary jo{x) and ji{x) (or ip{x) and ^/'(x)) the expression 
(I17p . in general, does not satisfy the boundary condition 

CI- 

Now we consider the model within the scope of which 
(as it has been discussed in Introduction) all the quan- 



tities given at the phase interface are represented by the 
linear combinations of the exponential functions. In par- 
ticular, let (p{x) and ip{x) be given by the expansions 

ip{x) ^Ao + Aie-^i"/2 ^ ^26-^^"/' + • • • , (20) 
ij{x) = Bo + Sie-Ti"/2 ^ ^ . . . ^ (21) 

where constants 7„ ^ 0, An and Bn will be specified in 
what follows. After the substitution of Eqs. (I^Dl) and (pij) 
in Eq. ([T7]) and the calculation of the integrals (details 
see in Appendix B), we obtain 



Jn{xt) 



n>0 



-lr.x/2 



exp 



2(l-a2)_ 



„(-) 



exp 



X 



2(1 



(22) 



and the notations have been introduced 

7i^^ = 7n + « ± + 2a7n + ^ ; (23) 

(24) 



2 



n I 



1 -F q;7„ 



52 



(25) 



(26) 



As it is seen from Eqs. ([20| and (f2T|) and Bo deter- 
minate the behavior of and ?/> (or jo and ji ) far from 



the system surface. Let us determinate the rest of the 
parameters 7„, An and Bn {n ^ 1) in such a way as 
to satisfy the balance condition (fT2|l and the boundary 
condition at the free surface ([13]). 



A. The determination of the parameters 



Now consider the boundary condition (fT3l) . Taking 

into account that 70 = 0, ^o = Q^j ^0 — ci-y 7o^^ — a ± a 
and using Eq. ([22]) . we have for an arbitrary i < 



J(x,t)|a 



Ao „ \ 2aH 
T-^°r"P2(T3^ 



Si 

T~^^^J^"P2(1 



7^^at 



2 t^i 



A2 
2 



<52 

B2—) exp 
1^2^ 2(1 



7r^«t 



exp 



/A2 , „ (^zX 



7^ ^at 
2(1^ 

72 ''ai 
2(1^ 



(27) 



If all the powers of the exponentials are different then 
J{0,t) = can be only at An = Bn = 0. However if 
each exponential function will appear in Eq. (|27p at least 
twice then this can lead to nonzero An and Bn- Bearing 
in mind this circumstance we determinate jn so that the 
equalities are held 



In fn— 1 



1,2,3,, 



(28) 



(+) 



in which 7 
be known [3 



(and respectively 7n-i) are considered to 
. Taking into account the notation 



and resolving the Eq 

„(+) 



in relation to 7„, one obtains 



(7«)l2 = 



At n = 1 and 7q 



(+) 



(l-a2) 

2a Eq. ([29]) gives 
4a 



(29) 



71 



1 



a^ 



The second value 71 = is the extraneous root of the 
Eq. ((29|) at n = 1. After the determination of 71 the 
values 7^^'' appearing in Eq. (j22p can be found from 
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TABLE I. The parameters 7„ appearing in Eqs. (f20)l and ([21}. 



n 





1 


2 


3 


4 


7n 

^(+) 

Jn 
Jn 




2a 



4q 


4a(3 + a2) 


8a(3 + a2)(l + a2) 


8a(l + a2)(a'* + 10a2 + 5) 


1-q2 
8a 


(1-q2)2 
2a{3 + a^f 


(l-a2)3 
32a(l +a2)2 


(l-a2)4 
2a(a* + lOa^ + 5)^ 


1 - q2 
2q 


(1-q2)2 

8a 


(1 - a2)3 


(l-a2)4 
32a(l + 0^)2 


l-a2 


(l-a2)2 


(l-a2)3 



Eq. Along similar a line one can obtain the val- 

ues 7„, 7^^^ for n > 1. In Table [J these values are given 
for n < 4. As it is seen from the table 7n (1 — a^)~", 
7^+' - (1 - a^)-" 7^^^ - (1 - The case of an 

arbitrary n is easily proved by induction using Eq. ([29]). 
Then considering the inequality 1 — ^ 1 we neglect 
by the exponentially small terms in the sum (|22p and re- 
strict ourselves by the terms with n ^ 3 only. To do this 
it will suffice to put 

4"^-^-S3— -0, A, = B„ = 0, n^4, (30) 

that corresponds to the first four terms of the expansions 
^ and (EH). 

Now let us define the rest of the nonzero constants An 
and Bn so that the condition ([?7|) holds. As a result, 



the constants i?„ are completely eliminated and we have 
finally 

J{x, t) = Ao(I - 6-^1-/2) exp 

1 — 

+{A^ + A3)(e-^^-/2 - e-^^-/2) exp tt^^^ 

(f — a'^y 

-KA(e-W2_,-W2)exp^ii±4^, (31) 

(1 — a'^y 

The simple but rather cumbersome calculations show 
that the flux defined by Eq. ((3T|) satisfies the initial con- 
ditions dni) and (fTe]) . 

It should be noted that the expression (PT|) can be rep- 
resented in the alternative form explicitly demonstrating 
the presence of the wave component in the mechanism of 
the solute transport 



J(x,t) 



aX 




aX 
AaX 



2 (l-a2)2 

72X a(3 + a2)2x 
(l-a2)3 



exp 



2 



4aA: 



exp 



722; 

2 



(l-a2)2 

a(3 + a2)2x 
(l-a2)3 



(32) 



where 



B. The solute concentration 



X = x + ai < 0; X 



at ^ 0: 



> 0,t ^ 0) 



As it is seen from Eq. (|32p the terms containing X can 
be considered as the concentration waves reflected from 
the surface and propagating to the interface. 
Except for the condition 

Ao + Ai+A2 + A:i= jo(0) - J(0, 0) = , 



following from Eqs. (|T8|) and ([20|) . the constants An ap- 
pearing in Eqs. (PT|) and ([5^ . remain up to now arbitrary 
and must be defined from other conditions that we shall 
consider in the next section. 



The solute concentration Cl in the liquid phase can be 
found in the same way as the flux has been defined. Let 



dCUxt) 
dt 



\t=~x/a = CLi(x) , 



and at the interface the relationships of the type of the 
expansions (PO)) and (pij) take place 

^e(x) = + af\-''^-/^ + af^e-''--'^ + • • • , 

where ipc(x) and ipc(x) are related to cl(x) and cli{x) by 
the equalities analogous to the Eqs. (fT5|) and (fT9)) . Then 
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FIG. 1. The concentration profiles for diflterent values of a — V/Vd and the surface solute concentration ci — Cs(0) in the 
solid. The material constants are Vd = 25(m/s), Vdi ~ 20(m/s), fee = 0.1, Co = 0.05. The solid lines are the solute profiles in 
the solid phase, the dash lines are the solute profiles in the liquid one. a) a = 0.95, ci = 0.06. In the inset the concentration 
profiles are shown when the interface is far from the surface. In both phases the curves are depicted by the solid lines. 
h) a — 0.85, ci = 0.06. In the inset the solute fiux in the liquid phase is shown. 



the solution of Eq. ([7]) (the latter is written down in the 
dimensionless form) has the form similar to the solution 
(j32p where A„ and _B„ must be replaced by and 6„. 
For determination of this parameters the mass conserva- 
tion law can be used. Substituting the expression (pS)) 
and the corresponding expression for the solute concen- 



tration Cl{x, t) into Eq. ([6]) and equating the coefficients 
at the linear independent functions one can write down 

and 6„ in terms An- As a result, using the found 
values for 7„ one has within the approximation as for 
Eq. ^ 



CL{x,t) = Co 

a 



exp 



2a 



(^2 + ^3)< exp 



aX 
l-a2 



exp 



aX 



(1 + 3a2)A3 
a(3-|-a2) 



exp 



2 (l-a2)2 
72X a(3-fa2)2x 



exp 



2 



4aX 



(l-a2)3 



exp 



(1 -a2)2 
722; a{i + a'^fX 



(l-a2)3 



where cq is the initial solute concentration in the melt 



Co — hm Cl{x, t) . 

t — > — oo 



At last, taking into account that at the interface 

cl(x) = Cl{x, t)\x=0 = ipcix) , 
one obtains for the coefficients of the expansion ([33] 



r 



Co 



a 



a 



(L) 
1 

(L) 
2 

(i) 



1 + 

2a 
l + a"' 



2a 

l-|-3a2 



{A2 + A3) - 



a{3 + a^) 



{A2 + A3) 
As- 



Ao 



l + 3a^ 
a{3 + a^) 



A. 



(34) 



(35) 



Now we consider the solute concentration in the solid 
phase. In line with the above we take the solute concen- 
tration cs{x) — Cs\x=o a-t the interface as 

(36) 
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where the second equahty means the condition of absence 
of the solute flux from surface of the solid usual for the 
Pick's diffusion. 

As regards the inhomogeneous partition coefficient 
ka(x), allowing for the Eqs. ([33]) and pG]) it can be pre- 
sented as 



fcoo(a) + fs{x,a)/c' 



1 + fL{x,a)/cl 



where 



(i = L,S). 



FIG. 2. The concentration profiles for a — 0.95, ci = 0.09, 
Co = 0.05, Vd = 25(m/s), Vdi = 20(m/s), fee = 0.1. The 
solid lines are the solute profiles in the solid phase, the dash 
lines are the solute profiles in the liquid one. 



Taking into account that 

J{x, t)\x=o =Ao- {Ao + A2+ ^3)6-^^^/2 

one substitutes the expression for cs{x), cl{x) and 
J(x,t)\x=Q in the interface boundary condition (|12p . 
Then using the equalities ([36)) . p3|) and ((35)l one obtains 
for the coefficients a? 



"0 



Cs = CO , 

2yln 1 - 



= -—-^r-{A2 + A,) 



AS) _ 



2a 
4(1 + a2) 
a(3 + a2) 



2a 

(^2 + ^3) 



(38) 



2(1 -«^) 
a(3 + a2)' 



A. 



A 



3 • 



As it is seen from Eqs. (|32j) and ([Ml)-® all the quanti- 
ties of interest to us are defined in terms of the param- 
eters Aq, A2 and A3. The value of the flux J£° = Aq 
at the interface infinitely removed from surface can be 
found from boundary condition (jl2p taken at cc ^ 00 



jr/a = cf-cT=co{l-l/koo), 

where ^00(0;) = jc^ is the partition coefficient for the 
infinite system given by ([2]) 



(\^a^)k, + a{yolVDi) 
fcoo(a) = \ ^ l-a'^ +aiyD/VDi) 



a < 1 
a ^ 1 



The constants A2 and A3 can be defined from the bound- 
ary conditions 



Cs(0)=ci 



dCs 



dx 



x=0 



0, 



III. DISCUSSION 

Figs. [iHS] present the solute distributions for cq — 0.05 
and different values of a = V/Vd and ci for some mo- 
ments of time t < 0. The solid lines are solute distri- 
butions in the solid phase, the dash lines are the solute 
distributions in the liquid. The jump of the concentra- 
tion takes place at the point of occurring of the interface. 
The terminal solute distributions Cs{x) at t = in the 
solid are given in Fig. [3l The behavior of the partition 
coefficient for different parameters is shown in Fig. 3] 

As it is seen from Figs.[T]-[3]the concentration curves are 
very sensitive to the variation of the parameters a and ci . 
In Fig.[T^ the concentration curves for a — 0.95 and ci — 
0.06 are shown. When the interface is far enough from 
the surface the solute distribution has the form typical 
for the stationary conditions of a solidification (inset in 
Fig. la). The motion of the interface near the surface 
gives rise to the increase of the solute concentration in 
the liquid phase. The curve at t = defines the terminal 
solute distribution in the solid (see also Fig. [3]). 

The partition coefficient for a = 0.95 and ci = 0.06 
is shown in Fig. and exhibits the monotonic increase 
reaching unity at the interface. 

At a = 0.85 and ci = 0.06 the solute distribution in 
the solid behaves not monotonically, having the maxi- 
mum removed from the surface (Fig. [T|d). In this case 
the partition coefficient exceeds unity in the near surface 
region (Fig. [4^). The reason may be in the following. 

At the high growth velocity V ^ Vd {a ^ 1) the 
interface reaches the surface practically simultaneously 
with the concentration disturbance moving with the ve- 
locity Vd ■ In this case the effects of the reflection have 
no time to affect markedly the value of the solute concen- 
tration at the interface. When the interface moves not 
enough rapidly (for example, at a = 0.85) the concentra- 
tion wave have time to reffect at the surface and to reach 
the solid phase in the near surface region. In this case 
the resulting solute flirx at the liquid side of the interface 
will be defined by the sum of two contributions: the flux 
of the solute atoms rejected by the interface and the op- 
positely directed flux of the the atoms reflected by the 
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FIG. 3. The terminal solute concentration profiles in the solid; co — 0.05, Vd ~ 25(m/s), Vdi = 20(m/s), fee = 0.1. 
a) ci = 0.06; b) ci = 0.09; a = V/Vd- 
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FIG. 4. The partition coefficient; Co = 0.05, Vd = 25(m/s), Vdi = 20(m/s), ke = 0.1. a) ci = 0.06; b) ci = 0.09; a = 



surface. In the end competition between them can lead 
to a decrease of the resulting flux and at a later time to 
its complete disappearance. 

In the inset in Fig.[T)D the solute flux in the liquid phase 
is shown. At the interface occurring at the sufflcient dis- 
tance from the surface the flux is directed to the liquid 
phase (t = -0.14,2:: « 0.12, J < 0). As the interface 
advances the flux at the interface practically disappears 
(at t « 0.1). The flux produced by the interface and the 
reflected flux will be equal in value. The jump of the 
concentration at this moment is absent and the partition 
coefficient reaches unity (Fig. 2^). At a later time owing 
to the reflected particle the resulting flux at the interface 
proves to be pointed toward the growth phase playing 
the role of the " external force" increasing the migration 
of the solute atoms through the interface to the solid (the 
inset in Fig. lb, t — —0.09, J > 0). As a result, this gives 
rise to rapid enough growth of the solute concentration 



at the solid side of the interface and the decrease at the 
liquid side of one (Fig. [T]d, t= - 0.04, J > 0). Beginning 
from the moment of time t w —0.1, the partition coeffl- 
cient becomes greater then unity (Fig.|4^), reaching unity 
only at the surface. 

At Fig. [3] the terminal solute distributions in the solid 
at different values a and ci are given. The concentration 
maximums are shifted to the surface with increasing ci 
(Fig. [Sb, ci = 0.09). At the same time at a = 0.95 
the partition coefficient does not behave monotonically, 
reaching the minimum near the surface. The following 
circumstance can be responsible for such behavior. 

When the interface moves with the velocity V ^ Yd 
the concentration disturbance has no time to propagate 
to considerable distances from the source of the distur- 
bance. As a result, the deviations from the initial con- 
centration Co will be only exhibited within the thin liq- 
uid layer near the surface. The greater is the velocity, 
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the thinner is this layer. The large enough value of Ci 
(relative to cq) presumes the drastic growth of the solute 
concentration in such a thin layer in comparison with 
the solid side of the interface. In Fig. [2] it is seen that 
when the interface moves in the near surface region the 
concentration jump increases leading to the decrease of 
the partition coefficient (Fig. At the smaller values 
of ci (ci = 0.06, a = 0.95) the interface advance is ac- 
companied by the decrease of the jump (Fig. la) and the 
monotonic behavior of (Fig. 2^). The same monotonic 
behavior of is appeared at ci = 0.09 but the smaller 
a. 



IV. CONCLUSION 

In the given work the one dimensional model of the 
solute diffusion during rapid solidification of the binary 
alloy in the semi-infinite volume is presented. Within the 
scope of the model it is supposed that the planar interface 
moves with a constant velocity, local equilibrium near 
the interface in the bulk of the liquid phase is absent and 
the diffusion flux is related to the concentration gradient 
by the Maxwell-Cattaneo equation. The latter gives rise 
to the telegrapher equation for both the flux and the 
solute concentration. The solution of this equation in 
the region occupied by the melt has been found by the 
Riemann method [33| . The basic assumption concerns 
to the behavior of the flux and the solute concentrations 
at the moving phase interface. Based on the model form 
of the partition coefficient given in the work [l^ it has 
been suggested that these quantities should be given in 
the form of the expansions ^ and ^ Such 
choice has made it possible to obtain in an explicit form 
the expressions for the flux and the solute concentrations 



in the both phases and to analyze different regimes of 
rapid solidification in the near surface region. 

Within the scope of the model it has been shown that 
different types of the concentration profiles can be real- 
ized depending on the external parameters (the growth 
velocity V and the solute concentration at the surface of 
the solid ci). As regards the partition coefficient, it may 
behavior not monotonically reaching both the maximum 
and the minimum in the near surface region. 

It should be noted that if the expansions ([^0]) . (PT|) 
contain the infinite sum of terms, then the expressions 
for the fiux ([32]) and the solute concentrations (p4)) and 
(1551) valid only for a — V/Vd close to unity. The equal- 
ities pop defining the number of the terms in the sum 
(l20l) and ([2T|) ( similar equalities hold for the expansion 
([55]) as well) have been written by recognizing that at 
a — > 1 the terms with n ^ 4 give the exponential small 
contributions to Eqs. dSH), §M and dM]), (01]). However, 
one can consider the models of solidification when the 
sums (j20|) and (j2ip contain a finite number of terms. In 
particular, if one defines 1^9(2;) from the outset in the form 

ip{x) =Ao + Aie-'^i^/2 A^e-^'"''!'^ + Aae"''^^/^ (39) 

and similarly for '^(x) and other quantities of interest, 
then when the first condition from the equalities ([50)) are 
held the expressions found for the fiux (|3ip . ([5^ and 
the solute concentrations ([34]) . ([36]) will give the exact 
solution for this model valid for any a < 1 including the 
case of Q — > that corresponds to local equilibrium and 
the Fick's diffusion. In the latter case it is easy to verify 
that the expressions ([32]) and ([34]) satisfy the diffusion 
equations. 

The terminal solute distributions in the solid for the 
model (1591) are given in Fig. [S^ for different values a and 
ci = 0.09. Comparing Fig. [5^ with Fig. [3)3 it is seen that 
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FIG. 6. The figures to the Riemann method. 



the solute distributions at the small a differ rather dras- 
tically from that at the large a. At the slow interface 
motion the concentration profiles does not behave mono- 
tonically at the given ci. The concentration maximums 
are several times bigger than the solute concentration at 
the solid surface. The behavior of the partition coefficient 
is presented in Fig. [Sb. The non-monotonic behavior is 
also exhibited. In the near surface region reaches its 
maximum exceeding unity. 

Another considerable assumption used in the work is 
that the interface remains planar during the whole of 
the solidification process. The problem of the stabil- 
ity of the planar interface has been attracting attention 
over several decades [U H^, [36-38]). It has been es- 
tablished that a planar interface is stable for both small 
enough growth velocities and high enough ones. More 
precisely there are the critical velocities Vc and Va so 
that if V < Vc and V > Va the interface is morpholog- 
ical stable against small perturbations of its form. In 
the region Vc < V < Va the planar interface is unstable 
and degenerates in different crystal patterns of cellular 
or dendritic morphologies. In the work [2l| it has been 
shown that for the isothermal solidification in the local 
nonequilibrium system Va < Vd ■ 

It has been noted however that the theoretical study 
of the stability problem is only based on using the steady 
state regime of the interface motion in a unbounded 
medium. The influence of the system surface on the 
stability of the planar front requires further investigation. 
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Appendix A: The Riemann method 

Let it be required to find the solution of the linear 
hyperbolic equation 



at2 



dJ_ 

at 



(A.l) 



satisfying the initial conditions given at the curve T : 
t = g{x) (see Fig. [6^) 

J\t=g(x) = joix) 

dj 

t=g{x) 



dt 



The substitution J 



-t/2 



u makes it possible to lead 



Eq. (|A.1|) to a more simple form 
d'^u d^u 1 

with the initial conditions 



U\t=g{x) 

du 



dt 



t=g(x) 



(A.2) 

(A.3) 
(A.4) 



The characteristics of Eq. (jA.2[) are the straight lines x ± 
t — const. According to the Riemann method 33] if the 
characteristics go through the point M and intersect with 
the curve F at the points P and Q , then the solution of 
Eq. (|A.2p at the point M can be represented as 



The authors thank P. Galenko for discussion of the 
work. 



u{M) 



PQ 



f dv ^ dv ^ 
u\ ^ — dti + -TT-dxi 
\OXi oti 



(A.5) 



The integral in (|A.5|) is taken along the curve F from P points P and Q. The Riemann function v(M,Mi) for 
up Q and up and uq are the values of u, taken at the 
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Eq. (IA.2|) has the form 

viM,Ah) = Jo(5v/(a^-a:i)2-(t-ti)2) , (A.6) 

where Jo{x) is the Bessel function of zero order and 
du/dx is calculated along the curve as 



du 
dx 



t=g(x) 



Lp[{x) - ipi{x)g'(x) . 



(A.7) 



The Riemann method for arbitrary linear hyperbolic 
equations one can find, for example, in (33| . 

Now consider the solution of Eq. (|A.2p in the region 
X > 0, t < 0, X = X + at < 0, when the initial data are 
given at the straight line t = —x/a (see Fig. [6^). Instead 
of (|A.3|) and (|A.4[) we have 



u\t=-x/a ^ 3oix)e "^Z^" 



du 
'dt 



t— — X I OL 



(^Jo+Ji)e-"/'" 



(A.8) 
(A.9) 



If the point M has coordinates (a;, t), so it is easy to show 
that the points P and Q have the abscissas respectively 
equal to 



xp = -- 



a(x + 1) 
1 — a 



Ci(x — t) 
1 + a 



(A.IO) 



Consider the integral term in Eq. (jA.5[) . Using 
Eqs. (jA.7[) - (|A.10[) and the fact that along the pathway 
of integration dti = —dxi/a, one has 



a{x — t) 
1 + a 



J y \aaxi oil/ ■> ti=-xi/a 



(A.ll) 



a{x + 1) 
1-a 



where the notations are introduced 

fix) ^jo{x), 

= hoix) - -j'aix) - ^ ji{x) . 
la a^ 

Furthermore using the Riemann function (|A.6[) . it can 
show that 



/■I dv dv \ 
\adxi dtiJ 



ti— — xi I a 



X 



J! 



l(^y{x-x,)^-{t + x,/a)^) 



(A.12) 



2q! y^[x - - {t + 
At last, after of the substitution the integral (jA.lip 
I 



into Eq. (jA.5|) and using the equality J — e^*/^u, one 
obtains the solution of the starting equation (|A.1[) . with 
added conditions (jlSp and in the form represented 
by Eq. ^ in Sec. HH 



Appendix B: The calculation of the integrals 

Substituting Eq. and (HIl) into we have 

J(x,t) = ^ J„(x,t), (B.l) 



where 



Jn{x,t) = -B„4i) + A„42) + ^{exp 



ain{x + t)+X 



2(1 -a) 



exp 



(XJnix -t) + X 



2(1 + a) 



(B.2) 
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and 



a{x — t) 
1 + a 



X 



a{x + 1) 
1-a 
a{x — t) 
1 + a 



't/2 



dx,e~-^/^'' 



a(x + 1) 
1 — a 



Sn = 



1 + ajn 



4(^^^ix-x,)^^it + x,/a)^) 
^ix~x,r-{t + x,/a)^ 



(B.3) 



(B.4) 



(B.5) 



The calculation J^^ 

Making the substitution in tlic integral (|B.3p 

2aX a{x + t) 

= — Vxi , 

1 — 1 — a 

we have (for convenience the index n is omitted) 



aX 
1 — a- 



■ exp 



X' 



2(1 -a) 



3, 



where the notations are introduced 

1 



d = j e-''Vo(/3V2(l-z))dz, 



X'^X+{^-l){x + t), 



(B.6) 



(B.7) 



aX 



<0, p 



X 



Vl - a" 



>0. (B.9) 



Consider the integral 3- Using the definition of the Bessel 
function 



Jn 



z — Z' 



2 - 



E 



(-l)"(/3/2)2'"(z^z2)r 

m! r(m + 1) 



where T(x) is the Euler gamma-function, one represents 
the integral l|B.7p in the form 



(-l)"(/3/2)2" 

— ' to! r(TO + 1) 



e-f'^iz- z^rdz (B.IO) 



Calculating the latter integral [3J|, one obtains 



1/2 



n=0 



(B.ll) 

where /^(x) is the modified Bessel function of the first 
kind. Furthermore, we use the equality [35| 



m=0 



t 

m 

|z| - \2t\ > 
In our case 



(B.12) 



|z|-|2t| = 



S\X\ 



2a(l-a2) 

(B-8) and instead of Eq. (|RTT|) we have 



{a^Y + 2a7 + a'') > 



(B.13) 



where 



■ y/ 0^72 + 2a7 + 



(B.14) 



Substituting the expression (jRTsI into Eq. (|R6l) and 
taking into account that Ii/2ix) = (2/7ra;)^/^ sinh(x), we 
obtain 



t/1 ^ — 



■ exp 



X' 



2(1 -a) 



exp 



a(l - i^n)X 



2^„(l-a2) 
I 



exp 



a(l + i/„)X 



25„(l-a2 



(B.15) 



where Vn = i^iSn). At last, substituting X' from and 
Eq. (IB.SP into Eq. (jB.lSP one has 



.... = |.-,,./^{«p 



r 7^+^^ 1 




r ik-'x 1 


— exp 


1.2(1 -a2)J 


L2(l-a2)J 



(B.16) 
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The calculated J, 



(2) 



we have (the index n is omitted) 



Consider the integral JA . After substitution of the 



variable in Eq. (|B.4p 



X 



X — Xi 



(B.17) 



j(2)^_Aexp 
4a 



X' 



aX 



2(1 -a) 26(1 ~a^) 
aX 



2 



i/2S_ 



, / 1 /I - // aX \2 



V ^ 1 — a 



-e 



aX 



// aX \2 



(B.18) 



V V 1 — a 



-e 



To calculate the integral (jB.lSp we consider the equality (jB.15[) , having previously made the substitution (jB.17|) into 
Jn"' . After reducing common factors, we have 



aX 
l-a2 



1 // aX \2 



2 V a2 V VI - a' 



-eO^-^sinh 



aA 



2(5(1 -a2 



1 -a2 

Differentiating the latter equation with respect to A, one obtains 

aX / 



4a 



1 - a2 J; 



I 2 V „2 



aA 
1 -a2 



1-a^ // aA \2 



a^ V V 1 — a 



cosh 



-e 

aX 



cosh ■ 



az^A 



25(1 -a2) 2(5(1 -a2) 



One multiplies the latter equality by 



exp 



A' 



aA 



2(1 -a) 2(5(1 -a2) 



and using Eqs. (|B.18|) . (|B.5|) and (|B.8|) . one has 

" 2 I P 2(1 

1 r ra7„(a; + ^) + A 



(l-a2) ^ "2(l-a2). 



l-a2)J 



2(1 - a) 



exp 



a"fn{x - <) + A 



2(1 + a) 



At last, substitute Eqs. (|B.16|) and (jB.2ip into Eq. (|B.2p and as a result we have 



(B.19) 



(B.20) 



(B.21) 
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2(1 -a2) 



exp 



2(1 -a2) 



(B.22) 



where 



^" - 2 

I 
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